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Comment on "Solution of Classical Stochastic 
One-Dimensional Many-Body Systems" 

In a recent Letter, Bares and Mobilia ^ proposed the 
method to find solutions of the stochastic evolution oper- 
ator H — Hq + -Ji?! with a non-trivial quartic term Hi. 
They claim, "Because of the conservation of probability, 
an analog of the Wick theorem applies and all multipoint 
correlation functions can be computed." Using the Wick 
theorem, they expressed the density correlation functions 
as solutions of a closed set of integro-differential equa- 
tions. 

In this Comment, however, we show that applicability 
of Wick theorem is restricted to the case 7 = only. To 
discuss this point, let us consider the generating function 
of correlation functions ^[^;t] = (x| exp[^^ ^qOg]e~-^*|p) 

introduced by Santos et al. Here aq(aj) is a q- 
mode annihilation (creation) operator and is a Grass- 
mann number with anticommutation relations {^g, a^] — 
{iqi Op} = {£,q,S,p\ = for all momentum indices p, q. (x| 
and \p) are the left vacuum and the initial state, respec- 
tively, belonging to the even sector of Hilbert space . 
For the initial state considered in Ref. ||] , the generating 
function Z is 


Z[^;0] = exp 


E 

.<3>0 


/i^ COt(f) 
1+^2 cot2(f) 


(1) 


with ji = p/{l — p) 1^. As mentioned in Ref. the 
Wick theorem at i = follows from the Gaussian form of 
the generating function. The applicability of Wick the- 
orem for later time can be tested by checking that the 


generating function is a Gaussian for later time or not. 

Now let us suppose that Z be a Gaussian with the 
form Z[^;io] = exp[X;g>o ^0)^9^-9] at time t = to. 
For simplicity we assume that / is an odd function 
of q as in Eq.(^. After an infinitesimal time incre- 
ment dt, the generating function becomes to + dt] = 


Z[^M + dt{x\[H, 


\p), where [,] denotes 


the usual commutation relation and we use the fact 
that {x\H = 0. Using the property of the left vacuum 
((x|aj — cot(|)(x|a_g, see Ref. Q) and the commutation 

relations [aj , Q] = -^^C, C] = 0, (C = expE, ^gOq]), 
we obtain the time evolution of the generating function. 
In the absence of the quartic term [i.e. 7 = 0), the gen- 
erating function evolves as 


Z[^; to + dt] = exp^ ^ ^,^-5 


X [J{qM) - dt{u[q) + D{q))f{qM) + '^edtsmq]\ , (2) 


where v{q) — uj{q) + 2esingcot(|) [0. This implies that 
Z remains a Gaussian at later time, if Z[(,; 0] is a Gaus- 
sian (see Eq.(|^)). As a byproduct, we obtain a differential 
equation (DE) for /. Since the DE for / preserves the 
parity in q, f is an odd function for all t. In the presence 
of the quartic term {i.e. 7 7^ 0), however, the situation 
should be modified significantly. After performing the 
commutation relation between quartic term and C, we 
find 


Z[i;to]-\x\[HiX]e-"'°\p) = 2^(cosgcos 


P 


l)cot(|)C,^_ 


'J{P, to)f{q, to) + 2 ^ sin (7 sinp cot(|)<^g$_,/(g, to)/(p, to) 


E sinp/(p, to) 


^ C0s((7i -t- q3)6{qi +q2 + qz + q4)f{q3,to)f{q4,to)^qi^q2£.q3£.q 


qi ,'J2,?3,94 


+ E [C0s(g2 +93) -C0s((7i +q3)]S{qi +52 + 93 + 94) COt(y )/((7i, to)/(q3, io)/(g4, to)Cgi^g2C936 


91, 92,93, 94 
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Although to] is a Gaussian initially, the quartic terms 
in Eq.(||) prevent it from remaining a Gaussian at later 
time, thus break the Wick theorem. 

If we ignore the quartic term in the generating func- 
tion, we indeed find the integro-differential equation of 
g{q,t){= /(g,t)cot(|)) as given by Eq.(7) in Ref. 0. 
Therefore, it may be considered as a Hartree-Fock type 
approximation ||^. Since the quartic term is irrelevant 
in the model considered, Eq. (7) in Ref. ^ predicts the 
correct (leading) long time behavior of the density and 
the correlation function. However, for the systems with 
the relevant non-quadratic terms in the evolution oper- 
ator, the validity of the method proposed in Ref. is 
doubtful. 
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